To appear in The Astrophysical Journal, 763:45, 2013 Jan 20 



1 



Chemical and Physical Conditions in Molecular Cloud Core 
DC 000.4-19.5 (SL42) in Corona Australis* 

E. Hardegree-Ullman\ J. Harju^'^, M. Juvela^, O. Sipila^, D. C. B. Whittet^, and S. Hotzel^-^ 

Abstract 

Chemical reactions in starless molecular clouds are heavily dependent on 
interactions between gas phase material and solid phase dust and ices. We have 
observed the abundance and distribution of molecular gases in the cold, starless 
core DC 000.4—19.5 (SL42) in Corona Australis using data from the Swedish 
ESO Submillimctcr Telescope. We present column density maps determined from 
measurements of C^^O (J = 2—1,1 — 0) and N2H+ (J = 1 — 0) emission 
features. Herschel data of the same region allow a direct comparison to the dust 
component of the cloud core and provide evidence for gas phase depletion of CO at 
the highest extinctions. The dust color temperature in the core calculated from 
Herschel maps ranges from roughly 10.7 to 14.0 K. This range agrees with the 
previous determinations from Infrared Space Observatory and Planck observations. 
The column density profile of the core can be fitted with a Plummer-like density 
distribution approaching n(r) ^ at large distances. The core structure deviates 
clearly from a critical Bonnor-Ebert sphere. Instead, the core appears to be 
gravitationally bound and to lack thermal and turbulent support against the pressure 
of the surrounding low-density material: it may therefore be in the process of slow 
contraction. We test two chemical models and find that a steady-state depletion 
model agrees with the observed C^®0 column density profile and the observed 
A''(C-^*0) versus Ay relationship. 
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1. Introduction 



Nascent solar systems inherit their chemical inventories from precursor molecular cloud cores. 
Within these cores, molecules containing the abundant CHON group of elements are of greatest 
importance to astrochemistry, particularly in the context of the possible emergence of habitable 
planets and life. CO is a dominant species within this group, because of both its abundance 



and the vital role it plays in chemical pathways leading to complex organic molecules (Tielens 
fc Hageii||1982 see Whittet et al.|[20TT and references therein for additional discussion). 



Gravitationally bound starless cores, referred to as "prestellar" cores (di Francesco et al. 



2007 Ward-Thompson et al. 2007), represent sites of imminent future star and solar-system 



formation. In this paper, we investigate the interplay of molecular gas and dust in the starless 
core DC 000.4—19.5 associated with the Corona Australis (CrA) molecular cloud complex. 
DC 000.4—19.5 has been identified in many surveys of the region and possesses numerous 
alternative designations: Cloud 42/S42/SLDN42 ( [Sandqvist fc Lindroos|1976[ ), Cloud B ( [Rossano 



1978), Condensation C (Andreazza & Vilas-Boas 1996), Core 5 (Yonekura et al. 1999), and 



Condensation CoA7 ( Vilas-Boas et al.]|2000 ), for example. Henceforth, the core will be referred 
to as SL42. SL42 is an isolated molecular cloud located at R.A. 19^10™16!3, Decl. -37°08'37", 
J2000.0, at a distance of 130 pc (Neuhauser & Forbrich 2008). No embedded young stellar 
objects (YSOs) have been detected within 5' of SL42, corresponding to the approximate radius 
of its core; Ha 16, the nearest YSO (Marraco & Rydgren 1981), is a weak-line T Tauri star 
( Batalha et al.|[l998 Gregorio- Hetem fc Hetem||2002[ ) 7^4 from the center of SL42. It is possible 
that the molecular cloud core which gave rise to Ha 16 was originally part of SL42 as evidenced 
by the present-day extension of the SL42 envelope toward Ha 16. The quiescent nature of SL42 
restricts chemical reactions within the cloud to low temperature gas phase molecular interactions 
and surface reactions on dust grains. 

Section [2] presents the data collection and reduction techniques for both radio and infrared 
observations of SL42. Section [s] describes the production of C^*0, N2H+, and H2 maps, the 
determination of visual extinction within SL42, and the estimates of the cloud's mass and 
dynamical state. Section [4] addresses the attempts to match chemical models to the observed 
properties of SL42, and Section [5] discusses the implications of our analysis. Finally, Section [6] 
summarizes our findings. 



2. Data 
2.1. SEST Observations 

Radio data were collected using the 15 m Swedish ESO Submillimeter Telescope (SEST) in 2003 
March. The molecular lines observed were C^^O J = 1 - (109.8 GHz), C^^^O J = 2 - 1 
(219.6 GHz), and N2H+ J = 1 - (93.2 GHz). A region of about 12' x 8', centered on the 
ISOSS source 19102-3708 (R.A. 19hl0™16!3, Decl. -37°08'37", J2000.0) was mapped in the two 
C^^O transitions simultaneously with the 3 and 1.3 mm superconductor-insulator-superconductor 
receivers, using a grid spacing of 40". A region of ~ 3' x 3' around the C^^O emission peak 
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was thereafter mapped in N2H+ J = 1 — and C^^O J = 2 — 1 using a 20" spacing. Figure [ij 
shows the coverage of the data field over a negative infrared Digitized Sky Survey (DSS^ image 
of SL42. 

All observations were performed in the frequency switching mode using an integration time 
of one minute per position. The two mixers used at the same time were connected to a 2000 
channel acousto-optical spectrometer (AOS) which was split in two bands of 43 MHz each. 
The AOS channel width corresponds to 0.14 km s"^ at 93 GHz, 0.12 km a"^ at 110 GHz, and 
0.06 km s~i at 220 GHz. The system temperatures were around 240, 260, and 450 K at 93, 110, 
and 220 GHz, respectively. The rms noise level attained was typically ^ 0.1 K at 3 mm and 
~ 0.2 K at 1.3 mm. 

The half-power beam width (HPBW) of the antenna was 47" at 110 GHz and 25" at 220 GHz. 
The pointing and focus were checked at 3—4 hr intervals toward circumstellar SiO maser line 
sources, and we estimate that the pointing accuracy was typically about 3". 

Calibration was done by the chopper wheel method. To convert the observed antenna 
temperatures (T^) to the radiation temperatures (Tr) the former were divided by the assumed 
source-beam coupling efficiencies (rjc) for which we adopted the main beam efficiencies {timb) 
of the telescope interpolated to the frequencies used, i.e., 0.73, 0.71, and 0.61 at 93, 110, and 
220 GHz, respectively. 



2.2. SEST Data Reduction 

The spectra were reduced using the Continuum and Line Analysis Single-dish Software (CLASS) 
program of the Grenoble Image and Line Data Analysis Software (GILDAS) package for the 
Institut do Radioastronomic Millimetrique (IRAM^ (See Figure [2] for example spectra.) The 
reduction involved the subtraction of the baseline using a polynomial fit and fitting Gaussian 
profiles to the detected lines to calculate the peak antenna temperatures (T^), local standard 
of rest (LSR) velocities (wlsr), and the full- width half-max (FWHM) line widths (Aw). The 
intensities were converted to the main-beam brightness temperature scale by Tmb = T^/vmb, 
where r]MB is the main-beam efficiency. Assuming uniform beam filling, Tmb represents the 
brightness temperature Tg. 

The N2H+ J = 1 — lines contain seven hyperfine (hf) components. To fit the hf structure 



we used the frequencies and the relative line strengths given in Caselli et al. (1995). Besides 
vlsr and Av, the hf fit also gives the total optical depth (rtot) of the line which is the sum of 
peak optical thicknesses of the seven components. The fit can be used to calculate the optical 
depth profile (r„) of the line (as a function of radial velocity v). The C^*0 lines do not have hf 
structure, and the lines could be well fitted with a single Gaussian. The lines are likely to be 
optically thin in SL42, a relatively isolated, small cloud. 



http://archive.stsci.edu/dss/acknowledging.html ©Anglo- Australian Observatory (AAO) /Royal Ob- 



ser vatory, Edinburgh (ROE) 



http : //www . iram . f r/IRAMFR/GILDAS 
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The C^^O and N2H+ column densities were determined from the reduced spectra by 
assuming that the rotational levels are populated according to the local thermodynamic 
equilibrium (LTE). For C^^O it was furthermore assumed that the lines are optically thin. 

When the excitation temperature (Te^) and the integrated optical depth (J T^dv) of the 
observed transition u I arc known, the column density of the molecules in the upper level 
{Nu) can be calculated from the formula 

Stt 1 f 
Nu = -^—F{Tui,Te^) I T^dv, (1) 

where A^i is the Einstein coefficient for spontaneous emission, A = c/f„; is the wavelength of the 
transition, T„; = hvui/ks, and the function F{Tui,T) is defined by 

For an optically thin line, as in the case of C^^O, the integral / Tydv can be estimated from 
the integrated intensity using the antenna equation: 

j TB{v)dv « Tui[F{Tui,T,,) - F{T^i,Tkg)] J r^dv. (3) 

For the derivation of the C^^O excitation temperature, the C^*0J=:2 — 1 spectral map 
was convolved to the same angular resolution as the J = 1 — map, and the T^.^ was calculated 
from integrated intensity ratio of the two lines using the antenna equation. 



I- 



The integrated r of the N2H+ line was obtained from hf fit results assuming a Gaussian 
profile: 

/tt 

Tydv = — ^==Aurtot- (4) 

The of N2H+ was calculated by substituting the Tg and r values at the line peak into the 
antenna equation: 

Tsiv) = Tui[F{Tui,T,^) - F{Tuung)]{l - e"""). (5) 

The total column densities {N^o\) including all rotational levels were estimated according to 
the LTE assumption: 

A^tot = ^e^"/'=^===Z(Te.), (6) 
9u 

where Z{Tex) is the rotational partition function at the temperature T^x, and S„ is the energy 
of the level u with respect to the rotational ground level. 



2.3. Herschel Data 

The 15' X 15' Herschel maps of SL42 used in this study were extracted from the extensive 
mapping of the CrA region collected by the Spectral and Photometric Imaging Receiver (SPIRE) 
and Photodetector Array Camera and Spectrometer (PACS) instruments as part of the Herschel 
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Gould Belt Survey ( Andre et al.||2010 see also http://www.herschel.fr/cea/gouldbelt/en/). An 

overview of the Herschel Space Observatory is given in Pilbratt et al. (20101. The SPIRE and 



PACS instruments are described in Griffin et al. (2010) and Poglitsch et al. (2010), respectively. 



The pipeline-reduced data from the SPIRE+PACS parallel mode survey of CrA are publicly 
available at the Herschel Science Archive]^ The original CrA maps cover a region of about 5° x 2° 
at the wavelengths 500, 350, and 250 fim (SPIRE), and 160 and 70 fj,m (PACS). The FWHM 
values of the Gaussians fitted to the beam profiles are, in order of decreasing wavelength, 37", 
25", 18", 12"xl6", and 6"xl2", the PACS beams being clearly elongated in observations carried 
out in the parallel mode. 

In order to calculate the distributions of the dust temperature (T^) and the optical depth (ta) 
of the dust emission, the Herschel maps at 500, 350, 250, and 160 /im were first convolved with a 
Gaussian beam to a resolution of 40" (FWHM) , and the intensity distributions were fitted pixel 
by pixel with a modified blackbody fmiction, I^^ « B^{Td)Tiy cx B^[Td)v^ , which characterizes 
optically thin thermal dust emission at far-IR and submillimeter wavelengths. The observed 
values of j3, the spectral index, are typically close to 2.0 (e.g., |Boulanger et al.|1996 Schnee et al. 



2010 



2003 



Juvela et al. 



Desert et al. 



2008 



2011 ), but higher values have been reported for cold regions (e.g., Dupac et al. 



Planck Collaboration et al. 2011 ). The recent study of the Taurus- Auriga 



molecular cloud found an average value of ~ 1.8 but again significant anticorrelation between 
dust temperature and the spectral index. In our analysis we fixed the dust opacity spectral index 
to a constant value oi ji — 2.0, which is used in several Herschel studies so the Td and A^(H2) 
maps presented in this work are directly comparable with previous results. This is consistent 
with the Planck Early Cold Core (ECC) estimate for SL42, /3 = 2.4 ± 0.6 (see Section [O]), 



especially considering that the Planck value is the estimate for the cold dust component rather 
than for the total dust emission from this source. The pipeline-calibrated in-beam flux densities 
were color corrected to monochromatic flux densities at the standard wavelengths in the course 
of Td fitting using the adopted shape of the source spectrum (a modified blackbody spectrum 
with (3 = 2.0) and the spectral response functions of SPIRE and PACS photometers for an 
extended sourcj^ In addition, the SPIRE pipeline in-beam flux densities were converted from 
point source calibration to extended source calibration before the temperature fitting and color 
correction. 

The resulting intensity map (/25o^im) at A=250 yitm and the Td map are shown in Figure |3] 
The optical depth map (r25o^m), calculated from T25o^m = hbOfim/ B250nm{Td), is proportional 
to the molecular hydrogen column density (Af(H2)). The column density can be obtained from 



7V(H2) 



'*250/jm"lH2 



(7) 



where k is the dust "opacity" , or absorption cross-section per unit mass of gas, and (= 2.8 
amu = 4.65 x 10"^''' g) is the average mass of gas per H2 molecule (assuming 10% He). The 
H2 column density errors were estimated by a Monte Carlo method using the la error maps 
provided for the three SPIRE bands and a 7% uncertainty of the absolute calibration for all four 



' http: //herschel . esac . esa. liit/Scl6nce_Archive . shtml 



'^See the SPIRE and PACS Observer's Manuals at 'http: //herschel . esac . esa. int/Documentation. shtml 
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bands according to the information given in SPIRE and PACS manuals. Different realizations 
of the iV(H2) map were calculated by combining the four intensity maps with the corresponding 
error maps, assuming that the error in each pixel is normally distributed. The iV(H2) error in 
each pixel was obtained from the standard deviation of one thousand realizations. For the dust 
opacity at 250 /im, we used the value /ilso^m — 0.085 cm^ g~^ which is a factor of 1.17 lower 



than the value k: 



250/jm 



0.1 cm^ g ^ consistent with the parameter C250 {— l/'«250Mm) given in 



Table 1 of Hildebrand (1983). This choice was motivated by the comparison between H2 colmnn 
densities from Herschel and visual extinctions from 2MASS as described in Section [3.21 



2.4. ISOPHOT and Planck Observations of Infrared Continuum Emission 

The Infrared Space Observatory photo-polarimeter (ISOPHOT) Serendipity Survey (ISOSS) 
recorded the 170 /im sky brightness when the satellite was slewing between two pointed 



observations. SL42 is one of the "cold cores" (Toth et al. 2000) detected in this survey. The 



dust temperature calculated by correlating the ISOPHOT 170 /xm and Infrared Astronomical 
Satellite (IRAS)/IRAS Sky Survey Atlas (ISSA) 100 /Ltm intensities is Td = ll.lto.s K. A detailed 



description of the method can be found in Hotzel et al. (2001 ). The ISOPHOT detector had four 



pixels with slightly different tracks. The width of the slew was about 3', and the effective angular 
resolution in the scan direction was 2'.2. The surface brightness distribution recorded by the two 
pixels which passed near the cloud center shows a Gaussian bump. It has an FWHM of 5f7 and a 
maximum intensity of 52 MJy sr""'^. Assuming spherical symmetry, the total 170 /im flux density 
of the core is 162 Jy. It sits on top of a 15 MJy sr~^ pedestal, which is likely to correspond to 
the low-density envelope around the core. Its position was used as the origin position for the 
SEST observations. 

SL42 can be found in the Planck Early Release Compact Source Catalogue (ERCSC|^as the 
object PLCKECC GOOO.37-19.51. The core is detected at the Planck frequencies from 100 to 
857 GHz (A = 3.0 - 0.35 mm). The coordinates of the 857 GHz maximum (PLCKERC857 
G000.36- 19.49) are at R.A. and Deck offsets (51", -9") from the SEST origin. The dust 
temperature and the emissivity index in the cold core are = 10.2 ± 0.6 K and /3 — 2.4 ± 0.6, 
respectively. These are calculated from the 857, 545, and 353 GHz fluxes, after the subtraction 
of the warm component (see the Explanatory Supplemenip^ . The ERCSC gives furthermore 
a total flux density /SgsyGHz = 271 ± 9 Jy and an angular size of 6'.9 x 6!l (FWHM) for the 
core. The Planck beam (FWHM) is 5.'1 x 3'.9 at 857 GHz. The beam size is similar at other 
submillimeter wavelengths. 



"See 
12 



http : //irsa . ipac . caltech . edu/applications/planck/ 



irsa. ipac . caltech. edu/data/Planck/release/ercsc_vl . 3/explaiiatory_supplemeiit_vl . 3 .pdf 
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3. Analysis 



3.1. C^^O and N2H+ Maps, Line Characteristics, and Column Densities 

Figure [4] shows the column density maps of the C^^O J = 2 — 1 and N2H"'' J = 1 — hnes. A 
maximum C^^O detection of 4.53 ± 0.09 (10^^ cm~^) occurs at positional offsets of (0",0"), and 
a maximum N2H+ detection of 13.40 ± 1.60 (lO^^ cm-^) occurs at positional offsets of (34", 25"). 
The C^^O J = 2 - 1 and N2H+ J = 1 - hue parameters and total N{C^^O) and iV(N2H+) 
toward both the mentioned maxima and the second brightest positions in each molecule are listed 
in Table [T] The spectra from these four positions are shown in Figure [2] 



3.2. Visual Extinction {Ay) 

Dust extincts light from background sources via absorption and scattering to make them appear 
redder and dimmer. We used two methods in order to estimate the amount of visual extinction, 
and thereby dust, in SL42. First, visual extinctions (Ay) were determined using Two Micron 



All Sky Survey (2MASS) JHK photometry as described in |Shenoy et aT] ( |2008| for field stars 
toward SL42 (marked in green in Figure [T]) . 

Figure [5] shows the color-color diagram for all of the field stars employed in this Ay survey. 
We included only those stars from the catalog without flags indicating contamination or poor 
photometric quality of the observations. To remove unreddened or anomalous stars from the 
sample, the following color constraints were applied: 

1.6{H -K) + 0.0 <{ J-H) < 1.6{H - K) + 0.53 (8) 

QA<{H^K). (9) 

The value of the reddening vector for the entire CrA cloud is 1.6, but the slope in the J — H 
versus H — K diagram does not change when only the SL42 region is considered. Stars outside 
of the parallel lines in Figure [5] were excluded since they could not be dereddened to the intrinsic 
color lines. 

It should be noted that the visual extinctions calculated from the 2MASS color-color diagram 
are conservative estimates. They were calculated using the minimum possible color excesses 
{E{J — K)) necessary to place the reddened field stars on the intrinsic color lines. Most of the 
reddening vectors of the field stars in our sample cross the intrinsic color lines in two places, 
once in a region of late spectral type and once in a region of early spectral type. It is likely 
that most of our field stars are of late spectral type, so the minimum extinction values we have 
assumed are correct. Even if this is not the case for all of the field stars, the maximum extinction 
values are higher at most by a factor of two, and accounting for a few early type stars should 
not significantly change the trend seen in Figure |6j detailed in Section [5j 

Figure [5] presents no clear evidence for previously unknown YSOs in SL42. This is consistent 
with the Herschel data, which confirm that SL42 is a starless core because the dust temperature 
decreases toward the center of the cloud. Ha 16 appears in Figure [5] as the red point below our 
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H — K color cutoff. Wilking et al. ( 1992 1 classify Ha 16 as a class II YSO by its spectral energy 
distribution (SED). 

The 2MASS method for estimating extinction is considered reliable because it measures its 
effects directly without making assumptions about the emissive properties of the intervening 
dust. However, this method is only useful in the outer regions of SL42 because the opacity 
in the core is so high that field stars become undetectable. In order to sample the cloud's 
extinction more completely, we used the canonical relationship N{ll2)/Av = 9.4 x lO^'^ cm~^ 



mag (Bohlin et al. 1978 Martin et al. 2012) as a second method to obtain dust extinction 



values using Herschel data. The calculation of A^(H2) from the Herschel intensity maps of SL42 
depends inversely on the assumed dust opacity (see Equation ([t])); therefore, it is important 
that the correct value of kappa is chosen. The 2MASS photometry provides an independent 
measurement of extinction at the edge of the cloud, so we adjusted ^fgQ^jj^ such that A^(H2) 
estimates from Herschel converted to Ay overlap the 2MASS data at low extinctions. This 
scahng gives a value of ^Igg^j^j = 0.085 cm^ g~^. 



3.3. Mass and Dynamical State 

The H2 column density map created from Herschel data as described in Section [23| can be used to 
estimate the cloud mass. By choosing the region above 7V(H2) = 3.3 x 10^^ cm~^ (corresponding 



to Ay = 3.5 when the relationship N{R2)/Av = 9.4 x 10^° cm^^ ^ag-^ of |BohUn et al] ( |1978[ ) 
is adopted) we obtain a mass of 21 Mq at the assumed distance of 130 pc. 

The derivation of the gravitational potential energy requires an assumption about the density 
distribution in the cloud. The H2 column density map has a sharp maximum at the offsets (46", 
11"), R.A. 19hl0'°20!2, Decl. -37°08'26", J2000.0, from the origin of our molecular line maps. 
We calculated the azimuthally averaged H2 column density profile using this position as the 
cloud center. The result of the calculation is shown in Figure [7j 

We also made an attempt to fit the azimuthally averaged column density distribution with 
a "Plummer-type" ( Whitworth fc Ward-Thompson||2001 Plummer 1911) outwardly decreasing 



number density profile using the following functional form: 



n(r) = when r < i?out, (10) 

1 + 4 



where n{r) is the number density of H2 as a function of radius r, no is the gas number in the 
cloud center, is the "fiat radius" representing the radial distance where the density increase 
toward the center levels off, and i?out is the radius inside which this density profile is supposed 
to be valid. We assumed that the cloud envelope outside this radius has a constant H2 column 
density that forms a pedestal on top of which we see the column density enhancement owing to 
the cloud core. 

The best fit to the data is obtained with the following parameter values: rj = 2.0, no = 
1.2 X 10^ cm~'^, ro — 15"0, and i?out = 300". The exponent value 77 = 2 means that the density 
distribution approaches n(r) oc r^^ for r ^ rQ. The average column density far from the cloud 
center, A^(H2) = 3.3 x 10^^ cm^^ (corresponding to Ay = 3.5) is used as the "pedestal" added to 
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N{Il2) obtained from the density model. The predicted column density profile (convolved with a 
Gaussian to 40" resolution to allow comparison with the observed one) is indicated as a dashed 
line in Figure [7] The core mass calculated from the adopted H2 density profile up to r = -Rout 
is Mcorc = 16 Mq. The 5 Mq difference between the model mass and the cloud mass from the 
N(H2) map is accounted for by the "pedestal" with N{H2) = 3.3 x 10^^ cm'^. 

The agreement between the observed (averaged) and predicted 7V(H2) distributions is 
reasonably good, and we used the model density distribution to calculate the gravitational 
potential energy of the core (i?grav) from the integral 

/■flout /■'■ 

i^grav = -IGn^Gml^ / rn{r)dr / y^n{y)dy ^ -1.0 x lO^^ J, (11) 
Jo Jo 

where G is the gravitational constant and {— 2.8 amu = 4.65 x 10^^* g) is the average mass 
of gas per H2 molecule (assuming 10% He, n{r) = riua (?"))■ In the kinetic energy estimate we 
assumed the cloud gas is isothermal with Tkin = 10 K and has a constant non-thermal velocity 
dispersion, ctn.t. = 260 m s~^. The latter value is calculated from the C^^O spectra in the core 
region. Using these values we obtain 

^kin = ^Meoro[c^ + CTn.T.] = 5.0 X lO^^ J, (12) 

where 

Cs = \/kBTy,in/fh (13) 

is the isothermal sound speed, m = 2.33 amu is the average particle mass of the gas. We see 
that according to the estimates above, the core satisfies the equation, iJgrav + 2i5kin — 0, which 
in the absence of external pressure would indicate virial equilibrium. However, our cloud model 
also includes a low density envelope which should be exerting pressure on the core, so it seems 
the cloud is lacking sufficient thermal and/or turbulent support against collapse 

The estimation of the kinetic and potential energies involves some uncertainties. The gas 
kinetic temperature is not directly measured, but we use the average dust temperature. Also, 
the H2 column densities from thermal dust emission are calculated using color temperatures and 
a constant dust opacity. The systematic errors related to these approximations are difficult to 



estimate. The color temperature is larger than the true temperature in the core center (Nielbock 



et al. 2012 Juvela & Ysard 2012 Malinen et al. 2011), and using the former temperatures we 



probably underestimate the H2 column densities and the mass. On the other hand, varying 
the temperature by 1 K and the velocity dispersion by one spectrometer channel changes the 
estimate of kinetic energy by 30 percent. Nevertheless, it is remarkable that our best effort 
estimates suggest a balance 2£'kin/|£'gi.av| ~ Ij showing that the core is gravitationally bound 
dMcKee &: Ostriker|[2007l). 



^^In the presence of an external pressure, the condition of virial equilibrium applying to a system within a 
surface S can be written as -Bgrav + 2_Eitin ~ SPextV = 0, where Pext is the pressure on the surface 5, and V is 
the volume of the system. A Bonnor-Ebert sphere always satisfies this condition. 
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4. Modeling 



We have attempted to reproduce the observed C^^O and N2H+ hne profiles (Figure [2| and the 
associated column densities (Table [l]) by calculating simulated radial abundance profiles for C^^O 
and N2H'^ using a chemical model and by using the abundance profiles as input for a radiative 
transfer program. We now discuss the modeling process and present the results. 



4.1. Chemical Modeling 

To simulate chemical evolution in the core, we divided the Plummer sphere discussed in Section 
3.3| into 40 concentric shells. Chemical evolution was then calculated separately in each shell, 
yielding radial abundance profiles as a function of time — the chemical model (and the general 



modeling procedure) used here is the same as the one discussed in detail in Sipila 



(2012). The 



gas phase chemical reactions are adopted from the OSU reaction file osu_ 03_200^ The initial 
gas phase abundances and the surface reaction set (i.e., activation energies for select reactions) 



are adopted from Semenov et al. (2010). For the chemical calculations, the visual extinction was 
calculated from the density profile as a function of core radius assuming Ay = 3 at the edge of 
the model core. In the calculations, it is assumed that all grains are 0.1 fim in radius and that 
the cosmic ionization rate is C = 1-3 x 10~^^ s~^. 

The upper left panel in Figure [s] shows the radial abundance profiles of C^*0 and N2H+ 
(with respect to total hydrogen density nn at three different time steps). Evidently the two 
species evolve very differently. C^^O is strongly depleted at the core center at later times but 
retains a fairly constant abundance at the core edge as the core grows older. N2H+ on the other 
hand increases greatly in abundance throughout the core as time progresses. The upper right 
panel of Figure [s] shows the ratio of the column densities of N2H+ and C^^O through the center 
of the model core as a function of time; the growth of the N2H+ abundance with respect to the 
C-'^^O abundance as the core grows older is especially prominent here. 

The C^^O column density map of the spherically symmetric model cloud has a ring-like 
maximum with a radius of ~ 40" (5000 AU). In the convolved (40" resolution) image, the 
slight depression at the core center has about 10% lower column density than the annular 
maximum. The corresponding N2II+ map shows a sharp peak toward the core center. The 
widths (FWHM) of the C^^O and N2H+ distributions are 200" and 60", respectively. The model 
can thus qualitatively explain the different extents of the C^^O and N2H+ emissions and the fact 
that these molecules can peak at different locations. 

The observed C^^O maximum does not lie, however, symmetrically around the A^(Il2) 
maximum but on its northwestern side. Possibly this can be understood in terms of the 
asymmetric H2 column density distribution which has a shoulder on the side where C^^O peaks, 
and decreases more rapidly in other directions (see Figure |4][ A)). 



*See 



http: / /www.physics. ohio-state.edu/~eric/ 
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4.2. Spectral Line Modeling 



We attempted to reproduce the observed spectra toward the N2H+ peak by calculating simulated 
N2H+ J = 1 — and C^^O J = 2 — 1 line emission profiles with a Monte Carlo radiative transfer 



program (Juvela 1997), using the chemical modeling results as input. The simulated spectra 
are highly dependent on the radial abundances and hence on the column densities of N2H+ and 
C-'^^O, which are both variable functions of time (Figure |8|. 

The column density ratio of C^^O and N2H+ toward the (34", 25") position (cf. Table [l]) 
corresponds roughly to i ^ 4 — 5 x Iff* years in the model. The lower left and lower right panels 
of Figure [s] show the simulated C^^O J = 2 — 1 and N2H+ J — 1 — line profiles, respectively, 
at i ~ 4.5 X 10^ years. The observed peak intensities (cf. Table [l]) are reproduced rather well 
by the model at this time step. However, the optical depths of the model lines are much smaller 
than the observed thicknesses. This is tied to the column densities predicted by the model; while 
the N2H+ and C^*0 column density ratio in the model corresponds to the observed ratio at 
t ~ 4.5 X 10^ years, the actual column densities for both species are at this time clearly smaller 
than observed. This suggests that the observed column densities correspond to a later stage in 
the model. 

To compare the observed and modeled column densities across the cloud we have plotted 
these in Figure |9] as functions of the visual extinction Ay, calculated from Herschel (using 
N{)l2)/ Ay — 9.4x 10^*^ cm^^ mag^^). The model distribution corresponds tot ^ 4.5x 10^ years. 
It is evident that the observed C^^O column density rises much more steeply with increasing Ay 
than predicted by the model. The model either undervalues the CO production through gas- 
phase reactions and desorption or uses too high estimates for the CO destruction rates (through 
photodissociation and accretion) in low-density gas. Alternatively, the discrepancy is caused by 
the fact that the static model adopted here does not take account of the increase of the density 
in case the core is contracting. We return to this issue in Section [5. 3| 



5. Discussion 



5.1. CO Depletion 



CO is expected to become depleted from the gas phase by adsorption onto grains on timescales 
less than the expected lifetimes of dense cores (e.g., Whittet et al.||2010 and references therein). 
In a detailed study of the Taurus dark-cloud complex, Whittet et al. (20101 compare gas phase 



and solid phase CO in lines of sight toward background field stars with known extinction values 
from 2MASS data. They find that (1) the total CO (gas -I- solid) column density agrees well 
with the general relation N = 2.0 x lQ^^{Ay — 2.0) cm~^ proposed by Kainulainen et al. (2006); 
and (2) depletion begins to become significant at extinctions Ay > 5 — 10 in the Taurus cloud. 



The plot of N{C^^O) versus Ay based on our observations and extinction estimates data 
from 2MASS and Herschel is shown in Figure |6] The straight line is the linear relation from 



Kainulainen et al. (2006), and the clear divergence of the data from this line at high extinctions 



{Ay > 15) is most likely caused by depletion. No data on ice-phase abundances are available 
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that would allow us to directly quantify the uptake of CO into grain mantles in SL42, but we 
note that the form of the variation of A^(C"'^^0) with Ay in Figure [g] is qualitatively very similar 
to that found by Whittet et al. ( 2010[ ) for Taurus. This variation can be described by a sigmoidal 
function of the form y = 02 + (ai — a2)/(l + {x/xoY). ^ of this form to our data for SL42 
is shown in Figure [6] and compared with the corresponding curve for Taurus from [Whittet et al.| 
(2010). We conclude that the iV(C^^O) versus Ay relation in SL42 follows the same general 



empirical law as that found for Taurus, but with lower levels of depletion. 

As previously mentioned, the Herschel Ay data in Figure |6] were scaled in order to match 
2MASS data at low extinctions. The scaling was uniform and made under the assumption that 



the dust opacity (K250;jm) remains constant throughout the cloud. Martin et al. ( 2012 ) found that 
this assumption is not necessarily true; in the dense molecular interstellar medium of the galactic 
plane, the dust opacity increases with decreasing temperature. Because our dust temperature 
maps are averaged along the line of sight, it is not possible to correct for the likely increasing 
dust opacity at the core of SL42. However, the dust opacity should vary at most by a factor of 
2 — 3 times from edge to core of SL42, and this would effectively compress Figure [6] along the 
X-axis, only at the highest extinctions. The general conclusion of CO depletion at the core of 
SL42 would be unaffected. 

The offsets in peak C^^O and N2H+ concentrations provide supporting evidence for CO 
depletion in the core of SL42. Offsets between C^^O and N2H+ are seen in many other starless 



cores, e.g., L1498, L1495, L1400K, L1517B, L1544 (Tafalla et al. 2002), and Barnard 68 (Lada 



et al.||2003[ ). The differentiation of N2H+ and CO originates from the main production pathway 



of N2H"'" in the reaction between N2 and Hj^. The N2 molecule is formed from relatively slow 
neutral-neutral reactions, and its abundance builds up later than that of CO. On the other hand, 
the H;^ ion increases strongly only after the depletion of CO in the densest parts of molecular 
clouds. The different behaviors of C-bearing and N-bearing species have been discussed, e.g.. 



by Flower et al. (20051 and by Fontani et al. (2012), see their Appendix A. At very late stages 



of chemical evolution N2H+ is also likely to become depleted, but previous observations suggest 
that even then the N2H+/CO abundance ratio remains high. One possible explanation suggested 



by Flower et al. (2005) is that a low sticking coefficient of atomic N could contribute to the 



availability of N2 also in highly depleted regions. 



5.2. Dynamics 



The Plummer-like density profile with three parameters, uq (the central density), tq (the 
characteristic radius of the inner region with a flat density distribution), and 77 (describing the 
density power-law far from the center) has been used in several previous studies for fitting the 



observed column density profiles of prestellar cores or starless filaments (Whitworth & Ward- 



Thompson 2001 Tafalla et al. 2002). Dust continuum maps have often suggested very steep 



(77 > 4) density gradients near the edges of prestellar cores (Whitworth & Ward-Thompson 



2001), but in the case of SL42 the azimuthally averaged column density profile can be fit with 



77 = 2 up to a radial distance where the core is melted into the envelope of low-column density 



gas. 
and L1400K. 



Tafalla et al. (20021 found similar outer profiles as in SL42 for the prestellar cores L1495 
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The average density profile obtained for SL42, n{r) — no/(l + (r/ro)^), can also be used to 
approximate a Bonnor-Ebert sphere (BES), which is a pressure-bound, isothermal, hydrostatic 
sphere. As evident from Figure [?) the structure of SL42 differs, however, clearly from a critical 



BES which has been frequently used to describe globules and dense cores (e.g., Alves et al. 



2001). The distribution, with a small flat radius tq compared to the outer radius of the core 
corresponds to a highly super-critical BES with ^out ~ 25 where the balance between forces is 
unstable against any increase of the outside pressure (dashed green line in Figure [?]). The overall 



agreement is, however, better with the adopted Plummer-like function with n(r) 
radii. 



at large 



Clear deviations from the critical BES structure have been previously seen in clouds with 



regular shapes other than SL42. Kandori et al. (2005) found in their NIR study of Bok globules 



that while several starless globules can be approximated by nearly critical BESs with ^out = 
6.5 ± 2, the column density profiles of star forming globules resemble those for super-critical 
BESs. The largest value of fout was found for the core L694-2 with in-fall features in molecular 
line spectra. The hydrodynamical models of Kandori et al. (2005) are consistent with these 



observations: collapsing cores mimic super-critical BESs with i^out increasing with time. 

A density distribution approaching the power-law n{r) cx at large values of r was found 
in the early numerical simulations of the collapse of an isothermal sphere by |Bodenheimer fc] 
Sweigart (1968) and Larson ( |1969 ), and the existence of a self-similar solution with this form 
was shown by Penston ( 1969| by analytic methods. This motivated Keto & Caselli (20101 to 
describe the collapse of a prestellar core by a sequence of BESs where the central density increases 
and the flat radius decreases while the outer radius is constant and the density profile in the 
outer parts keeps its scaling as (see their Figure 1). The density distribution approaches 
that of a singular isothermal sphere {n{r) cx r~'^ throughout) which is the initial configuration 



in the inside-out collapse model of Shu (1977). 



Keto & Caselli (20101 predict the flat radius tq to be equal to the product of the sound 



speed and the free-fall time in the core center (see their Equation (2)). Allowing for turbulent 
velocity dispersion we may write a slightly modified formula for rg: 



32Gpo 



(14) 



where po is the central density. Substituting the adopted values of the parameters, Tkin = 10 K, 



ctn.t. = 260 m s 



no 



1.2 X 10® cm we get for the expected flat radius tq — 1900 AU or 



15", i.e., the same value as determined from the observations. 

The calculated density profile and the agreement between the flat radius with the theoretical 
prediction suggest that SL42 is indeed contracting. According to the model of |Keto &: CaseUT] 
(2010) the contraction speed is highest around the flat radius (15" in our case) and clearly 
subsonic at this stage of the evolution. The detection of infall is not possible with the resolution 
and sensitivity of the present data. 
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5.3. Modeling 



Our pseudo-time dependent chemistry model discussed in Sections |4 . 1 1 and |4 . 2 1 failed to reproduce 
the C^^O column densities in the outer parts of SL42. Judging from the deduced density profile, 
the cloud is evolving dynamically, and probably the main difficulty here is to find realistic initial 
conditions for the chemistry taking the earlier history of the cloud into account. In their genuinely 



time-dependent model for a contracting core, Keto & Caselli (2010) found that CO reaches its 



equilibrium abundance in a short time compared with the dynamical timescale. It seems therefore 
justified to examine if a steady-state model can explain the observed C^^O column density profile. 

Assuming a steady state, the fractional C^^O abundance can be described by a simple 
function of the gas density. The accretion and desorption rates (cni^^ s^^) in an isothermal 
cloud can be written as ^n(H2)n(Ci*0, gas) = An(Il2)^X(Ci^0, gas) and Bn{C^^O,dust) = 
(112)^5^(0^^0, dust), respectively. The adsorption coefficient A and the desorption coefficient 



B depend on the dust properties, temperature, and the fiux of heating particles (e.g., Leger 



eTaT] [19851 [Hasegawa fc Herbst| [19931 [Caselh et all [20021 [Hotzel et ah] [20021 [Keto fc Caselli 
2008 2010 1. The coefficients A and B have the following relationships with the desorption 
and depletion timescales, tqn and tqff: defined by [Keto fc C^^e^ ( [2008| : An{R2) = I/ton, 
B — 1/toff- The definitions of tqn and tqff are given in Equations (8) and (11) of Keto & 



Caselh (2008). 



In what follows we denote the no-depletion fractional abundance of C^^O by Xq {— 
X(C^*0,gas) -I- X(C^^O, dust)), and the gas-phase fractional abundance simply by X {— 
X(C-^*0, gas)). The steady-state condition, An(R2)^X = Bn{ll2){Xn — X), imphes the relation 
( [Keto fc Caseiii[[2008[ Equation (12)): 



X 



1 + Toff /ton 



(15) 



Using the notation of Hotzel et al. (2002), Equation (8), this equation becomes X = Xo/{l + 
^n(H2)). In the idealized situation where the gas isothermal, the grain size distribution remains 
unchanged, and the same desorption mechanisms are effective in the whole volume considered, 
the ratio A/B is constant, and X is a function of n{ll2) only. 

In the steady-state model, several (rather uncertain) parameters affecting the gas-phase 
C^^O abundance degenerate into one quantity, the ratio A/B with the dimension cm~'^. This 
ratio can be therefore identified with the inverse of the characteristic density, say riD, at which 
depletion owing to freezing out becomes substantial. With this definition the equation for X can 
be rewritten as 

^0 



^ l + n(H2)/nD' 
where (using the notation of Keto fc Caselli[[2008 ) is 

3 X 10-7CHjlO-i7exp(-Eco/70) 



SoRde-o-Vr 



(16) 



(17) 



where is the cosmic ray ionization rate of H2, -Eco is the binding energy of CO onto ice, Sq 
is the sticking coefficient of CO to the dust in a collision, the product RdgCr is the grain surface 
area per H2 molecule, and Vr is the thermal speed of C^^O molecules. 
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Assuming an isothermal cloud at 10 K, and using the dust parameters and the cosmic ray 
ionization rate listed in Keto & Caselli (20081, no obtains the value ^ 1.2 x 10"* cm 

"^^ cm^ (per H2) and 



^. Here the 
^ s-i. The 



assumed total grain surface area is ~ 1.4 x 10~ cm (per H2) and = 3 x 10^ s 
density can be made higher by, e.g., increasing the minimum grain size (which decreases the 
total grain surface area and reduces accretion) or by increasing the cosmic ray ionization rate 
(which enhances desorption). 

We have plotted in Figure [9] the C^^O column density as a function of Ay as predicted by 
the steady-state model (red solid line), together with the iV(C^*0) versus Ay from observations 
(crosses) and the profile predicted by our time-dependent chemistry model (black dashed line) 
which were discussed in Section 4.1 The Ay axis of the model data (lines) is calculated from 



the spherically symmetric cloud model described in Section [3^31 For the steady-state depletion 
model a rough agreement with the observed column density profile was attained by setting 
Xq{C^^O) = 2.4 X 10"'', riD = 3.8 x 10^ cm'^. The change in ud with respect to the value quoted 
above corresponds to a decrease in the total grain surface area by a factor of three, or to an 
increase in the cosmic ray ionization rate (or the overall efficiency of desorption mechanisms) 
by the same amount. A decrease in the effective total grain surface area can be justified by 
assuming that CO is efficiently desorbed from the smallest grains by other processes than cosmic 
ray heating (Hotzel et al. 2002 Section[6]). We note that Keto & Caselli (2010) needed to increase 
the desorption rate by an order of magnitude from that implied by the direct cosmic ray heating 
(denominator of the formula for tt-d) to obtain agreement with the observed CO spectra toward 
L1544. 

The prediction of the steady-state model agrees better with the observed CO column density 
profile than the time-dependent model. We note that while the time-dependent model was 
optimized to match with the observed C^^0/N2H+ column density toward the core center, we 
have no estimate for the N2H+ column density profile based on the steady-state model, hence we 
cannot in this case attempt to simultaneously match C^^O and N2H+ profiles to observations. 
Thus, we cannot draw conclusions as to which model is better in the sense of reproducing all of 
the observed data. 

The agreement between C^^O observations and the steady-state model is not perfect as the 
model column density does not seem to rise steeply enough. Furthermore, the observations and 
the model predictions at low extinctions appear to have a horizontal offset of about 2 magnitudes 



from the empirical linear relationship of Kainulainen et al. ( 2006 ) . This offset can probably be 



accounted for by the uncertainty related to the zero levels of the pipeline reduced Herschel 
maps. The calibration can be improved in the future when the Planck data of the region become 
publicly available. 

The prediction of the steady-state model gives equally good agreement with the observed 
C^'^O column density profile as the sigmoidal function. This suggests the characteristic shape of 
the A^(C^^O) versus A^ correlation seen in SL42 and previously in Taurus reflects the density 
structure in the core envelopes. In the steady-state model, local variation in the CO depletion 
depends mainly on the density. 

Both the density structure and the CO depletion law found in SL42 agree with the model 



of Keto & Caselli (2010) for a thermally super-critical prestellar core. The outer density profile 
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follows the power law n(r) cx r~^, the flat radius is consistent with the estimated physical 
parameters in the core center, and the depletion time scale is short compared with the contraction 
time scale so that the CO abundance distribution does not deviate much from that for a static 
cloud with the same physical characteristics at any moment. 



6. Summary 

We have studied the structure and chemistry of the relatively isolated prestellar dense cloud 
Sandqvist & Lindroos 42 in Corona Australis using the dust emission maps observed with the 
Herschel satellite and C^'^^O and N2H+ line maps observed with the SEST radio telescope. The 
H2 column density map as calculated from Herschel data shows a sharp peak suggesting a high 
central density. The azimuthally averaged column density distribution can be fitted with the 
characteristic prestellar density profile which scales as in the outer part. 

C^^O and N2H+ peak at different locations, about 50" (6500 AU) apart. The N2H+ 
distribution is compact, and the maximum coincides with the A'^(H2) peak. The separation 
between the N2H+ and C^^O can be probably understood in terms of CO depletion. According 
to chemistry models, the N2H"'" abundance builds up slowly and is benefited by the disappearance 
of CO in the densest parts of molecular clouds. 

The relation between A^(C^^O) and Ay (or iV(H2)) is clearly not linear. The curved 



distribution of points on a A^(C^^O) versus Ay plot is similar to that found by Whittet et al. 



( |2010| ) for Taurus and can be fitted with a sigmoidal function. The curvature is most likely 



caused by accretion of CO onto grain surfaces. 

The shape of the iV(C^^O) versus Ay relation can be explained fairly well with a steady- 
state depletion model, where the fractional C^^O abundance is a simple function of density: 
A:(C^'*0) Xo(C^*0)/(l + n/no), where Xo(C^^O) is the no-depletion fractional abundance 



and riD is a characteristic density depending on the desorption and adsorption rates (Keto & 
Caselli|[2008l [Hotzel et al.|[2002| . 



The density structure and the CO depletion law observed in SL42 agree with the models of 



Keto & Caselli (2008 20101 for the evolution of a thermally super-critical prestellar core, where 
the subsonic contraction is described by a series of nearly hydrostatic cores with decreasing flat 
radii and density profiles in the outer parts. Our molecular line data are not appropriate for 
detecting the possible radial infiows. 
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Fig. 1. — Data coverage of the sample over a negative infrared DSS image of SL42. The red field is 
the distribution of C^*0 measurements, and the size of the circles corresponds to the HPBW of the 
antenna at 220 GHz (25"). The blue diamonds show the distribution of N2IT'' measurements. Locations 
where the observational noise exceeded the emission signals (i.e., S/N<1 or S<0) are not shown. The 
green squares indicate the locations of stars with good quality 2MASS JHK photometry which met the 
condition H -K > 0.4. The image center is located at R.A. 19''10"'16':3, Decl. -37°08'37", J2000.0. 
The image is 15' x 15'. 
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Fig. 2. — Example spectra from the sample. The left column shows C^^O J = 1 — and J = 2 — 1 

emission lines. The right column shows N2H+ J = 1 — hues. The four positions represent the brightest 
spots in C^**0 (top two rows) and N2H+ (bottom two rows). The multiple peaks in N2H+ spectra 
are due to hyperfine splitting. The origin position for the SEST observations, (0",0"), is located at 
R.A. 19''10'"16!3, Decl. -37°08'37", J2000.0. 
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Fig. 3. — Colorbar units are MJy sr^^, K, and cm^^, respectively. (A) Herschel intensity map 
of thermal dust emission at the wavelength A = 250 /im in a 15' x 15' region around SL42. 
(B) Dust temperature map of the same region, created by fitting a modified blackbody function 
to the Herschel intensity maps at 500, 350, 250, and 160 /im (see Section 2.3 1. Blue contours 
range from 11 K at the center to 15 K in increments of 1 K. (C) H2 column density map 
determined from the far-infrared optical depth (see Equation Q). The black contour shows 
where 7V(H2) = 3.3 x lO^^ cm'^. 
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Fig. 4. — Contour maps of (A) C^'^^O and (B) N2H+ superimposed on the H2 column density map 
of SL42. The contours for C^^O (red) range from 1 to 4 (10^^ cm~^) in increments of 1; yeUow 
indicates a depression contour. The contours for N2H+ (blue) range from 2 to 11 (10^^ cm^^) in 
increments of 3. The maximum C^^O detection is located at (0", 0"), and the maximum N2H+ 
detection is located at (34", 25") giving a spatial offset of ^ 42" (^ 5500 AU at a distance of 
130 pc). The black contour shows where -/V(H2) = 3.3 x 10^^ cm~^. The colorbars are in units 
of cm~^. 23 



Table 1: Line Characteristics and the C^^Q and N2H+ Column Densities in Selected Positions 







Ci^O J 


= 2-1 








(Aa, AS) 


A, peak 




Av 


T 


T 

ex 


N{C^^O) 


(") 


(K) 


(km s"-'^) 


(km s"'^) 




(K) 


(IQis cm-2) 


(0,0) 


3.52 (0.07) 


5.53 (0.01) 


0.72 (0.01) 




6.65 (0.28) 


4.53 (0.09) 


(-25,34) 


3.66 (0.07) 


5.45 (0.01) 


0.70 (0.01) 




6.77 (0.39) 


4.47 (0.11) 


(34,25) 


3.17 (0.06) 


5.48 (0.01) 


0.78 (0.01) 




6.87 (0.38) 


4.18 (0.09) 


(59,-9) 


3.01 (0.06) 


5.52 (0.01) 


0.77 (0.01) 




7.40 (0.41) 


3.51 (0.06) 






N2H+ J 


= 1-0 








(Aa,A<5) 


A, peak 


V-LSR 


A?; 


T 


T 

ex 


Af(N2H+) 


(") 


(K) 


(km s-i) 


(km s~^) 




(K) 


(IO12 cm-2) 


(0,0) 


0.63 (0.07) 


5.67 (0.02) 


0.52 (0.04) 


4.75 (1.98) 


4.14 (0.32) 


5.97 (2.62) 


(-25,34) 


0.51 (0.06) 


5.50 (0.02) 


0.73 (0.05) 




5a 


2.35 (0.45) 


(34,25) 


1.28 (0.05) 


5.66 (0.01) 


0.52 (0.01) 


8.00 (0.89) 


4.98 (0.11) 


13.40 (1.60) 


(59,-9) 


0.99 (0.06) 


5.70 (0.01) 


0.50 (0.02) 


9.03 (1.36) 


4.29 (0.11) 


11.60 (1.87) 



Note. " r and T^^ could not determined. The N2H+ column density was calculated by assuming 
optically thin emission with T^x = 5. 
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H-K (mag) 

Fig. 5. — Color-color diagram of all reddened {H — K > 0.4) field stars behind SL42 (equivalent to 
the green squares in Figure [l]) Green circles denote stars within the SEST data field, while yellow 
circles denote stars beyond the SEST field. Ha 16 is shown as a red point below the H ~ K cutoff. 
The reddening slope for SL42 is 1.6. Red and blue mark the intrinsic color lines for giants and dwarfs, 
respectively. 
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Fig. 6. — Top: C^*0 column density as a function of extinction. Red points: Interpolated A'^(C^*0) 
and field star extinctions from 2MASS JHK photometry. Green points: SEST N{C}^0) and extinctions 
from Herschel maps using N(R2)/Av = 9.4 x 10^°cm~^ mag~^ from Bohlin et al. (19781. Blue line: 
sigmoidal fit to data from the Taurus cloud (see Whittet et al. 2010 for details). Red line: average 
fit based on observations of molecular clouds with low depletion (A^ = 2.0 x 10^*(j4x/ — 2.0) cm~^) 
( Kainulainen et al.|2006"| . Black line: sigmoidal fit to SL42 data (iV = 4.07953+ (0.14037-4.07953)/(l + 
(^v/11. 27966)^ '^"^^) cm-2). Bottom: SEST Af(N2H+) and extinctions from Herschel maps. Blue 
points: A'^(N2H''') found by fitting the hyperfine structure of the spectral feature. Orange points: 
A'^(N2H"'") found assuming optically thin conditions. 
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Fig. 7. — Azimuthally averaged H2 column density distribution of SL42 as a function of 
radial distance in arcseconds from the maximum (crosses), which occurs at R.A. 19'^10™20?2, 
Decl. -37°08'26", J2000.0, or offsets from the SEST origin of (46", 11"). The solid orange curve 
shows the H2 column density calculated from a spherically symmetric density profile following 
the function n{r) — no/{l + (r/ro)^), where uq = 1.2 x 10® cm^'^, and tq = 15". The solid green 
curve represents a critically stable Bonnor-Ebert sphere (^out = 6.45) with the maximum column 
density agreeing with the observed one. The observed distribution can be fitted fairly well also 
with a highly super-critical (^out ~ 25) Bonnor-Ebert sphere (dashed green curve). 
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Fig. 8. — Upper left: C^^O (red) and N2H+ (green) radial abundances with respect to total 
hydrogen density nn in the chemical model at i = 4 x IC* years (solid lines), t — 8 x 10* years 
(dashed lines) and t = 2 x 10^ years (dotted lines). Upper right: the ratio of N2H+ and C^*0 
column densities through the center of the core model as a function of time. Lower left: the 
simulated C^^O J = 2 — 1 line emission profile at t = 4.5 x lO"' years in the antenna temperature 
scale. Lower right: the simulated N2H+ J = 1 — line emission profile at t = 4.5 x 10"' years in 
the antenna temperature scale. 
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Fig. 9. — Chemical modeling fits to SL42 observations. Green diamonds correspond to observed 
data. The dashed black line (top panel) and black crosses (bottom panel) show the best fits of 



the Sipila (2012) model, while the red line shows the best fit from the iKeto & Caselli (2008) 



model. The sigmoidal function fitted to the C^*0 data is shown in blue, and the empirical fit to 
undepleted clouds ( Kainulainen et aL]|2006 1 is in orange. 
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